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Abstract 

vo ; 

The fundamental nature of mass is one of the greatest mysteries in physics. The Higgs mechanism is a theoretically 
appealing way to account for the different masses of elementary particles and implies the existence of a new, yet 
unseen particle, the Higgs boson. We search for evidence of a light scalar (e.g. a Higgs boson) in the radiative decays 
of the narrow T(3S) resonance: T(3S) — > "fA°, A — > //+ . Such an object appears in extensions of the Standard 
Model, where a light CP-odd Higgs boson naturally couples strongly to 6-quarks. We find no evidence for such 
processes in a sample of 122 x 10 6 T(3S) decays collected by the BaBAR collaboration at the PEP-II B-factory, and 
set 90% C.L. upper limits on the branching fraction product B(T(35) -> 7A ) x B(A° n + n~) at (0.25-5.2) x 10" 6 
in the mass range 0.212 < m^o < 9.3 GeV. We also set a limit on the dimuon branching fraction of the r/b meson 
^ ! B(j] b — > [i + [>r) < 0.8% at 90% C.L. The results are preliminary. 
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I. INTRODUCTION 

The concept of mass is one of the most intuitive ideas in physics since it is present in everyday human experience. 
Yet the fundamental nature of mass remains one of the greatest mysteries in physics. The Higgs mechanism is a 
theoretically appealing way to account for the different masses of elementary particles [l|. The Higgs mechanism 
implies the existence of at least one new particle called the Higgs boson, which is the only Standard Model (SM) 
particle yet to be observed. If it is found, its discovery will have a profound effect on our fundamental understanding 
of matter. A single Standard Model Higgs boson is required to be heavy, with the mass constrained by direct searches 
to m H > 114.4GeV Q and m H ^ 170 GeV [4], and by precision electroweak measurements to m H = 129±4 4 GeV 0]. 

The Standard Model and the simplest electroweak symmetry breaking scenario suffer from quadratic divergences in 
the radiative corrections to the mass parameter of the Higgs potential. Several theories beyond the Standard Model 
that regulate these divergences have been proposed. Supersymmetry Q is one such model; however, in its simplest 
form (the Minimal Supersymmetric Standard Model, MSSM) questions of parameter fine-tuning and "naturalness" 
of the Higgs mass scale remain. 

Theoretical efforts to solve unattractive features of MSSM often result in models that introduce additional Higgs 
fields, with one of them naturally light. For instance, the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) 0] introduces a singlet Higgs field. A linear combination of this singlet state with a member of the 
electroweak doublet produces a CP-odd Higgs state A whose mass is not required to be large. Direct searches typi- 
cally constrain m^o to be below 2mb |8J] making it accessible to decays of T resonances. An ideal place to search for 
such CP-odd Higgs would be T — > 7A , as originally proposed by Wilczek Q. A study of the NMSSM parameter 
space predicts the branching fraction to this final state to be as high as 10 -4 . 

Other new physics models, motivated by astrophysical observations, predict similar light states. One recent exam- 
ple [III proposes a light axion-like pseudoscalar boson a decaying predominantly to leptons and predicts the branching 
fraction B(T — * 7 a) to be between 10 -6 -10 -5 [ll|. Empirical motivation for a light Higgs search comes from the 
HyperCP experiment [T^]. HyperCP observed three anomalous events in the £ — ► final state, that have been 

interpreted as a light scalar with mass of 214.3 MeV decaying into a pair of muons [13j |. The large datasets available 
at BABAR allow us to place stringent constraints on such models. 

If a light scalar A exists, the pattern of its decays would depend on its mass. Assuming no invisible (neutralino) 
decays [13], for low masses m^o < 2m T , relevant for the HyperCP [l2j and axion [ll| interpretations, the dominant 
decay mode should be A — > Significantly above the tau threshold, A — > t + t~ would dominate, and the 

hadronic decays may also be significant. This analysis searches for the radiative production of ^4° in T(3S) decays, 
with A decaying into muons: 

r(3S) -> j A ; A -> n+p- 

The current best limit on the branching fraction B(T — > jA ) with A° — > comes from a measurement by 

the CLEO collaboration on T(1S) [H|. The quoted limits on B(T(1S) -> 7A ) x B(A Q -> n+ >~) are in the range 
(1-20) xl0~ 6 for m^o < 3.6 GeV. There are currently no competitive measurements at the higher-mass T resonances 
or for the values of m^o above the tt threshold. 

In the following, we describe a search for a resonance in the dimuon invariant mass distribution for fully reconstructed 
final state Y(3S) — > 7(/i + ^~). We assume that the decay width of the resonance is negligibly small compared to 
experimental resolution, as expected [111, [l6| for m^o sufficiently far from the mass of the rjb We also assume 

that the resonance is a scalar (or pseudo-scalar) particle; while significance of any peak does not depend on this 
assumption, the signal efficiency and, therefore, the extracted branching fractions are computed for a spin-0 particle. 
In addition, following the recent discovery of the r]b meson in T(3S) decays [17J, we look for the leptonic decay of the 
775 through the chain T(3S) — + 777b, r)b — * If the recently discovered state is the conventional quark-antiquark 

rjb meson, its leptonic width is expected to be negligible. Thus, setting a limit on the dimuon branching fraction sheds 
some light on the nature of the recently discovered state. We assume T(rjt,) = 10 MeV, which is expected in most 
theoretical models and is consistent with BaBAR results [TtJ • 

II. THE BABAR DETECTOR AND DATASET 

We search for two-body transitions T(3S) — > 7A , followed by decay A — > in a sample of (121.8 ± 1.2) x 10 6 

Y(3S) decays collected with the BABAR detector at the PEP-II asymmetric-energy e + e~ collider at the Stanford Linear 
Accelerator Center. The data were collected at the nominal center-of-mass (CM) energy E cm = 10.355 GeV. The CM 
frame was boosted relative to the detector approximately along the detector's magnetic field axis by (3 Z = 0.469. 

We use a sample of 78.5 fb _1 accumulated on T(AS) resonance (T(AS) sample) for studies of the continuum back- 
grounds; since T(4S) is three orders of magnitude broader than T(3S), the branching fraction T(4S) — > 7A is 
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expected to be negligible. For characterization of the background events and selection optimization we also use a 
sample of 2.4 fb" 1 collected 30MeV below the T(35) resonance. 

Since the BABAR detector is described in detail elsewhere 18], only the components of the detector crucial to 
this analysis are summarized below. Charged particle tracking is provided by a five-layer double-sided silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH). Photons and neutral pions are identified and measured using the 
electromagnetic calorimeter (EMC), which comprises 6580 thallium-doped Csl crystals. These systems are mounted 
inside a 1.5-T solenoidal superconducting magnet. The Instrumented Flux Return (IFR) forms the return yoke of the 
superconducting coil, instrumented in the central barrel region with limited streamer tubes, and in the endcap regions 
with the resistive-plate chambers, for the identification of muons and the detection of clusters produced by neutral 
hadrons. We use the GEANT software to simulate interactions of particles traversing the BABAR detector, taking 
into account the varying detector conditions and beam backgrounds. 

III. EVENT SELECTION 

We select events with exactly two oppositely-charged tracks and a single energetic photon with a CM energy 
E* > 0.5 GeV. We allow other photons to be present in the event as long as their CM energies are below 0.5 GeV. We 
assign a muon mass hypothesis to the two tracks (henceforth referred to as muon candidates), and require that they 
form a geometric vertex with the x^tx < 20 (for 1 degree of freedom) , displaced transversely by at most 2 cm from the 
nominal location of the e + e _ interaction region. We perform a kinematic fit to the T(3S) candidate formed from the 
two muon candidates and the energetic photon, constraining the CM energy of the Y(3S) candidate, within the beam 
energy spread, to the total beam energy y/s. We also assume that the T(3S) candidate originates from the interaction 
region. The kinematic fit improves the invariant mass resolution of the muon pair. We place a requirement on the 
kinematic fit Xr(3S) < ^9 (for 6 degrees of freedom), which corresponds to the probability to reject good kinematic 
fits of less than 10~ 6 . The kinematic fit x 2 > together with a requirement that the total mass of the T(3S) candidate is 
within 2 GeV of y/s, suppresses background events with more than two muons and a photon in the final state, such as 
cascade decays T(3S) — > 7 Xb(2P) — * jjT(1S) — > etc. We further require that the momentum of the dimuon 

candidate A a and the photon direction are back-to-back in the CM frame to within 0.07 radians, and select events in 
which the cosine of the angle between the muon direction and A direction in the center of mass of A is less than 
0.88. We reject events in which neither muon candidate is positively identified in the IFR. 

The kinematic selection described above is highly efficient for signal events. After the selection, the backgrounds 
are dominated by two types of QED processes: "continuum" e + e~ — > 7/i + fi~ events in which a photon is emitted in 
the initial or final state, and the initial-state radiation (ISR) production of the vector mesons J/tp, il>(2S), and T(1S), 
which subsequently decay into muon pairs. In order to suppress contributions from ISR-produced p° — > ir + ir~ and 
<f> — > K + K~ final states in which a pion or a kaon is misidentified as a muon or decays {e.g. through K + — > /x+z^), 
we require that both muons are positively identified when we look for A candidates in the range m^o < 1.05 GeV. 
Finally, when selecting candidate events in the r/b mass region m w ~ 9.39 GeV, we require that no secondary 
photon above a CM energy of E% = 0.08 GeV is present in the event; this requirement suppresses decay chains 
T(35) -> j 2 Xb(2S) -> 7i7 2 r(LS'), in which the photon 72 has a typical CM energy of 100 MeV. 

We use Monte Carlo samples generated at 20 values of m^o over a broad range 0.212 < m^o < 9.5 GeV of possible 
A masses to measure selection efficiency for the signal events. The efficiency varies between 24-44%, depending on 
the dimuon invariant mass. 

IV. EXTRACTION OF SIGNAL YIELDS 

The invariant mass spectrum for the selected candidates in the T(3S) dataset is shown in Fig. [TJ We extract the 
yield of signal events as a function of the assumed mass m^o in the interval 0.212 < m^a < 9.3 GeV by performing a 
series of unbinned extended maximum likelihood fits to the distribution of the "reduced mass" 

m R = \J m 2 ^ - Ami ' W 

The choice of this variable is motivated by the distribution of the continuum background from e + e~ — > 7 /i + ^i~, 
which is a smooth function of mn across the entire range of interest, in particular, the region near the kinematic 
threshold ps 2m^ (tor ~ 0). Each fit is performed over a small range of around the value expected for a 
particular m^o. We use the T(AS) sample to determine the probability density functions (PDFs) for the continuum 
background in each fit window, which agree within statistical uncertainties with Monte Carlo simulations. We use 
a threshold (hyperbolic) function to describe the background below jur < 0.23 GeV; its parameters are fixed to the 
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FIG. 1: Distribution of the dimuon invariant mass ?ri„+„- in the T(3S) data. Black histogram shows the distribution for the 
selection in which only one of two muons is required to be positively identified. The peak from e + e~ — » 7isrP°(770), p° — » tt + iy~ , 
in which one of the pions is misidentified as a muon, is clearly visible. Blue (lower) histogram shows the distribution for the 
selection in which both muons are positively identified. The ISR-produced peaks at J/ip and T(1S) masses are visible. 



values determined from the fits to the T(4S) dataset. Elsewhere the background is well described in each limited ran 
range by a first-order (ma < 9.3 GeV) or second-order > 9.3 GeV) polynomial. 

The signal PDF is described by a sum of two Crystal Ball functions [2fJ] with tail parameters on either side of the 

maximum. The signal PDFs are centered around the expected values of ttlr = <J m 2 A0 — 4m^ and have the typical 

resolution of 2— lOMeV, which increases monotonically with m^o. We determine the PDF as a function of m A o using 
a set of high-statistics simulated samples of signal events, and we interpolate PDF parameters and signal efficiency 
values linearly between simulated points. We determine the uncertainty in the PDF parameters by comparing the 
distributions of the simulated and reconstructed e + e~ — > 7isrJ/0, J/0 ~ * events. 

Known resonances, such as J/0, 0(25*), and T(1S), are present in our sample in specific intervals of m^, and 
constitute peaking background. We include these contributions in the fit where appropriate, and describe the shape of 
the resonances using the same functional form as for the signal, a sum of two Crystal Ball functions, with parameters 
determined from the dedicated MC samples. We do not search for A signal in the immediate vicinity of J/0 and 
0(2Sj, ignoring the region of w ±40MeV around J/0 (approximately ±5cr) and ps ±25MeV (« ±3cr) around ip(2S). 

For each assumed value of m^o, we perform a likelihood fit to the mn distribution under the following conditions: 

• 0.212 < m^n < 0.5 GeV: we use a fixed interval 0.01 < iur < 0.55 GeV. The fits are done in 2 MeV steps in m A o. 
We use a threshold function to describe the combinatorial background PDF below ttir < 0.23 GeV, and constrain 
it to the shape determined from the large T(4S) dataset. For rriR > 0.23 GeV, we describe the background by a 
first-order Chebyshev polynomial and float its shape, while requiring continuity at nriR = 0.23 GeV. Signal and 
background yields are free parameters in the fit. 

• 0.5 < m A o < 1.05 GeV: we use sliding intervals /i — 0.2 < inn < M + 0.1 GeV (where /i is the mean of the 
signal distribution of m^). We perform fits in 3 MeV steps in m^o. First-order polynomial coefficient of the 
background PDF, signal and background yields are free parameters in the fit. 

• 1.05 < m^o < 2.9 GeV: we use sliding intervals fi — 0.2 < mji < /i + 0.1 GeV and perform fits in 5 MeV 
steps in m^o. First-order polynomial coefficient of the background PDF, signal and background yields are free 
parameters in the fit. 

• 2.9 < m A o < 3.055 GeV and 3.135 < m A a < 3.395 GeV: we use a fixed interval 2.7 < m R < 3.5 GeV; 5 MeV 
steps in m A o. First-order polynomial coefficient of the background PDF, signal, J/0, and background yields 
are free parameters in the fit. 
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• 3.4 < m^o < 3.55 GeV: we use sliding intervals /i — 0.2 < itir < fi + 0.1 GeV and perform fits in 5 MeV 
steps in m^o. First-order polynomial coefficient of the background PDF, signal and background yields are free 
parameters in the fit. 

• 3.55 < m^o < 3.66 GeV and 3.71 < m^a < 4.0 GeV: we use fixed interval 3.35 < run < 4.1 GeV; 5 MeV steps in 
m^o. First-order polynomial coefficient of the background PDF, signal, tp(2S), and background yields are free 
parameters in the fit. 

• 4.0 < m^o < 9.3 GeV: we use sliding intervals /i — 0.2 < itir < fi+ 0.1 GeV; 5 MeV steps in m^o. First-order 
polynomial coefficient of the background PDF, signal and background yields are free parameters in the fit. 

• 1% region (m Vb = 9.390 GeV): we use a fixed interval 9.2 < ran < 9.6 GeV. We constrain the contribution from 
e+e~ — > 7isrT(15') to the expectation from the T(4S) dataset (436 ± 50 events). Background PDF shape 
(second-order Chebyshev polynomial), yields of Y(3S) —> 7%b(2P) — > 'fjT^S), signal T(3S) — > 777;, events, and 
background yields are free parameters in the fit. 

The step sizes in each interval correspond approximately to the resolution in m^o . 

As a crosscheck, we also perform a set of two-dimensional maximum likelihood fits to the joint distribution of mn 
and cos6>* M , where is the CM polar angle of the dimuon pair. For scalar A°, the angular distribution of the signal 
events is expected to be 1 + cos 2 0*^, modulo acceptance and efficiency effects. The distribution of the most dominant 
QED backgrounds is more strongly peaked in the forward and backward directions. Thus, the angular distribution 
can help distinguish any signal peaks from the QED background. We find that the results of the two-dimensional fit 
are consistent with the one-dimensional fit to m# only. 

V. SYSTEMATIC UNCERTAINTIES 

The largest systematic uncertainty in B(T(3S) — > 7^4°) comes from the measurement of the selection efficiency 
We compare the overall selection efficiency between the data and the Monte Carlo simulation by measuring the 
absolute cross section da/dniR for the radiative QED process e + e~ — > 7/z + /z~ over the broad kinematic range 
< mji < 9.6 GeV, using a sample of 2.4fb _1 collected 30 MeV below the Y(3S). We use the ratio of measured 
to expected cross sections to correct the signal selection efficiency as a function of m^o. This correction reaches up 
to 20% at low values of m^o. We use half of the applied correction, or its statistical uncertainty of 2%, whichever 
is larger, as the systematic uncertainty on the signal efficiency. This uncertainty accounts for effects of selection 
efficiency, reconstruction efficiency (for both charged tracks and the photon), trigger efficiency, and the uncertainty 
in estimating the integrated luminosity. We find the largest difference between the data and Monte Carlo simulation 
in modeling of muon identification efficiency. 

We determine the uncertainty in the signal and peaking background PDFs by comparing the data and simulated 
distributions of e + e~ — * JisrJ/iP events. We correct for the observed 24% difference (5.3 MeV in the simulations 
versus 6.6 MeV in the data) in the width of the ran distribution for these events, and use half of the correction to 
estimate the systematic uncertainty on the signal yield. This is the dominant systematic uncertainty on the signal 
yield for ra^a > 0.4 GeV. Likewise, we find that changes in the tail parameters of the Crystal Ball PDF describing 
the J/ijj peak lead to variations in event yield of less than 1%. We use this estimate as a systematic error in the signal 
yield due to uncertainty in tail parameters. 

We find excellent agreement in the shape of the continuum background distributions for ijir < 0.23 GeV between 
T(3S) and T(4S) data. We determine the PDF in the fits to T(AS) data, and propagate their uncertainties to the 
T(3S) data, where these contributions do not exceed a(B) = 0.3 x 10~ 6 . For the higher masses mn > 0.23 GeV, the 
background PDF parameters are floated in the likelihood fit. 

We test for possible bias in the fitted value of the signal yield with a large ensemble of pseudo-experiments. For 
each experiment, we generate a sample of background events according to the number and the PDF observed in the 
data, and add a pre-determined number of signal events from fully-reconstructed signal Monte Carlo samples. The 
bias is consistent with zero for all values of ra^o , and we assign a branching fraction uncertainty of a(B) = 0.02 x 10~ 6 
at all values of m^o to cover the statistical variations in the results of the test. 

The uncertainties in PDF parameters of both signal and background and the bias uncertainty affect the signal 
yield (and therefore significance of any peak); signal efficiency uncertainty does not. The effect of the systematic 
uncertainties on the signal yield is generally small. The statistical and systematic uncertainties on the branching 
fraction B(T(3S) — ► 7^4°) as a function of ra^o are shown in Fig.O 
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FIG. 2: Statistical and systematic uncertainty on the product of branching fractions B(T(3S) — > 7^4°) X B(A° — * p + p~) as 
a function of m A a, extracted from the fits to the T(3S) data. Statistical errors are shown as red dot-dashed line, systematic 
uncertainties are shown as blue dotted line, and the total uncertainty, computed as a quadrature sum of statistical and systematic 
errors, is the solid black line. The shaded areas show the regions around the J/tp and ip(2S) resonances excluded from the 
search. 



VI. STATISTICAL INTERPRETATION 



In the event of a positive observation, statistical significance of a particular peak needs to be assessed. Convention- 
ally, this is done by computing the likelihood ratio variable 



S{m A o) = v /21og(L max /L ) (2) 

where -L max is the maximum likelihood value for a fit with a floated signal yield centered at m^o , and Lq is the value 
of the likelihood for the fixed zero signal yield. Under the null hypothesis (no signal events in the data), the signed 
quantity, sign(N s i g )S is expected to be normal-distributed (where iV s ; g is the fitted signal yield). The distribution for 
our T(35 l ) dataset is shown in Fig. [31 Since our scans have 0(2000) m^o points, we should expect several statistical 
fluctuations at the level of S w 3, even for a null signal hypothesis. 

For a single experiment consisting of Atrial uncorrelated measurements, the probability to observe S > 6> max and 
iV s ig > is 

P(5 max ; Atrial) « AftrialPGSmax; 1) = N t , lal Erfc (W^ (JJ) 

(this approximation is good for P(iS max ; Atrial) <C !)■ The "trial factor" Atrial is difficult to estimate analytically. 
Instead, we determine Atrial by inspecting results of 10 8 Monte Carlo pseudo-experiments; for each experiment we 
generate 1951 random values Xi according to 

Xi — Xi-!p + r\/l - p 2 (4) 

where p = 0.84 is the average correlation coefficient between adjacent "bins" Xi and Xi-i, as determined by T(3S) 
fits, and r is a normal-distributed random number. We then compute the maximum value <S max for each pseudo- 
experiment. The cumulative distribution of <S max describes the chance that a pure background would fluctuate to 
produce a signal peak at any value of m^o with the likelihood ratio variable S > 5 max . From the fit to the distribution 
of 5 max from 10 8 pseudo-experiments, we determine the effective trial factor Atrial = 1436 ± 16. Finally, we convert 
the likelihood ratio variable S to the true statistical significance S in terms of gaussian a: 

S' w x/^Erfc- 1 [Ar trlal ErfcGS max /\/2)l {S > 4) (5) 
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FIG. 3: Distribution of the likelihood ratio variable <S (with additive systematic uncertainties included for the fits to the Y(3S) 
dataset. The blue curve is the Gaussian fit with fixed fi = and a = 1. 

In particular, the threshold for observing "evidence" for signal at any value of m^o is S > 3.0; this corresponds to 
the likelihood ratio at a particular value of m^o of S > 4.8. 

VII. RESULTS AND CONCLUSIONS 

For a small number of fits in the scan over the T(3S) dataset, we observe local likelihood ratio values S of about 
3o\ The most significant peak is at = 4.940 ± 0.003 GeV (likelihood ratio value S = 3.0, including systematics; 
B = (1.9±0.7±0.1) x 10~ 6 ). The second most-significant peak is at m A o = 0.426±0.001 GeV (likelihood ratio value 
S = 2.9, including systematics; B = (3.1 ± 1.1 ± 0.3) x 10~ 6 ). The plots for these points are shown in Fig. [4] and 
Fig. El The peak at m^o = 4.940 GeV is theoretically disfavored (since it is significantly above the r threshold), while 
the peak at m^o = 0.426 GeV is in the range predicted by the axion model [llj. Neither of the peaks are significant, 
however, when we take into account the trial factor discussed in Section PVTl At least 80% of our pseudo-experiments 
contain a fluctuation with S = 3<r or more. 

Since we do not observe a significant excess of events above the background in the range 0.212 < m^o < 9.3 GeV, 
we set upper limits on the branching fraction B(T(3S) — ► 7^4°) x B(A° — > /i + /i~). We add statistical and systematic 
uncertainties (which include the additive errors on the signal yield and multiplicative uncertainties on the signal 
efficiency and the number of recorded T(3S) decays) in quadrature. The 90% C.L. Bayesian upper limits, computed 
with a uniform prior and assuming a Gaussian likelihood function, are shown in Fig. as a function of mass m^o . 
The limits fluctuate depending on the central value of the signal yield returned by a particular fit, and range from 
0.25 x 10~ 6 to 5.2 x 10~ 6 . 

We do not observe any significant signal at the HyperCP mass, m^o — 0.214 GeV (see Fig. [TU|) . We find B(T(3S) — > 
7 A°(214)) = (0.12t°oil ± 0.17) x 10- 6 , and set an upper limit B(T(3S) -> 7 v4°(214)) < 0.8 x lO" 6 at 90% C.L. 

A fit to the ri b region is shown in Fig. HQ We find B(T(3S) -> 7%) x B(rj b -> = (0.2 ± 3.0 ± 0.9) x 10~ 6 , 

consistent with zero. Taking into account the BABAR measurement of B(T(3S) — > 777b) = (4.8 ± 0.5 ± 1.2) x 10 -4 , 
we can derive B{rj b -► ^ + M") = (0-0 ± 0.6 ± 0.2)%, or an upper limit B^t -> fi + ^) < 0.8% at 90% C.L. This is 
consistent with expectations from the quark model. All results above are preliminary. 

The limits we set are more stringent than those reported by the CLEO collaboration recently [l5| • Our limits rule 
out much of the parameter space allowed by the light Higgs [T3| and axion [ll| models. 



13 




> 

o 

m 
o 



120 



100 



80 



60 



40 



20 



Preliminary 



X 2 /df= 105.8/118 




4.75 



4.8 



4.85 



4.9 



4.95 



m R (GeV) 



FIG. 4: The fit for m^o = 4.940 GeV in Y(3S) dataset. The bottom graph shows the run distribution (solid points), overlaid 
by the full PDF (solid blue line). Also shown are the contributions from the signal at m A o = 4.940 GeV (solid red line) and the 
continuum background (dashed black line). The top plot shows the normalized residuals p — (data — fit)/<r(data) with unit error 
bars. The signal peak corresponds to the likelihood ratio variable S = 3.0, including systematics, and B = (1.9±0.7±0.1) x 10~ 6 . 
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FIG. 5: The fit for m A o = 0.426 GeV in Y(3S) dataset. The bottom graph shows the run distribution (solid points), overlaid 
by the full PDF (solid blue line). Also shown are the contributions from the signal at m A o = 0.426 GeV (solid red line) and the 
continuum background (dashed black line). The top plot shows the normalized residuals p — (data — fit)/<r(data) with unit error 
bars. The signal peak corresponds to the likelihood ratio variable S = 2.9, including systematics, and B = (3.1±1.1±0.3) x 10 -6 . 
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FIG. 6: Upper limits on the product of branching fractions B(Y(3S) — > 7^4°) X B(A° — > ) as a function of m A o from the 

fits to T(3S) data. The shaded areas show the regions around the J/i/j and -0(25) resonances excluded from the search. 
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FIG. 7: Upper limits on the product of branching fractions B(T(3S) — > x B(A° — ► ) as a function of m A o in the 

range 0.212 < m A o < 1.05 GeV from the fits to T(3S) data. 
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FIG. 8: Upper limits on the product of branching fractions B(T(SS) —* 7 A ) x B(A° — * fi + fi~) as a function of rn A o in the 
range 1 < m A o < 4 GeV from the fits to T(3S) data. The shaded areas show the regions around the J/ip and ip(2S) resonances 
excluded from the search. 
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FIG. 9: Upper limits on the product of branching fractions B(T(3S) 
range 4 < m A o < 9.3 GeV from the fits to Y(3S) data. 
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FIG. 10: The fit for m A o = 0.214 GeV (HyperCP candidate) in the Y(3S) dataset. The bottom graph shows the ma distribution 
(solid points), overlaid by the full PDF (solid blue line). Also shown are the contributions from the signal at m A o — 0.214 GeV 
(solid red line) and the continuum background (dashed black line). The inset zooms in on the HyperCP region of interest. The 
top plot shows the normalized residuals p = (data — nt)/cr(data) with unit error bars. 
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FIG. 11: The fit for the r/b region in Y(3S) dataset. The bottom graph shows the van distribution (solid points), overlaid 
by the full PDF (solid blue line). Also shown are the contributions from the signal at at m Vb = 9.389 GeV (solid red line), 
background from the e + e~ — > 7isr2"(1£') (dot-dashed green line), background from T(35) — > 7Xt(2P), Xfc(2P) — » jT(1S) 
(dotted magenta line), and the continuum background (dashed black line). The inset shows the signal, e + e~ — > 7isrT(1S), 
and Xfc(2P) — > yT(lS) in more detail. The top plot shows the normalized residuals p = (data — fit)/cr(data) with unit error 
bars. 
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